Hydrodynamic excitations of trapped dipolar fermions 
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A single-component Fermi gas of polarized dipolar particles in a harmonic trap can undergo a mechanical 
collapse due to the attractive part of the dipole-dipole interaction. This phenomenon can be conveniently ma- 
nipulated by the shape of the external trapping potential. We investigate the signatures of the instability by 
studying the spectrum of low-lying collective excitations of the system in the hydrodynamic regime. To this 
end, we employ a time-dependent variational method as well as exact numerical solutions of the hydrodynamic 
equations of the system. 
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Following a spectacular series of experimental achieve- 
ments with Bose-Einstein condensates in dilute trapped 
atomic gases (for a review see, e.g., 01 01), dilute gaseous 
atomic Fermi systems have drawn intensified interest only 
very recently. After the first successful realization of quan- 
tum degeneracy in a potassium gas Q], several other groups 
have cooled fermionic 6 Li below the Fermi temperature and 
trapped it magnetically or optically |4|, |5j |6J . Very recently, 
two-species Fermi-Bose mixtures have been also achieved in 
6 Li- 23 Na H and in 40 K- 87 Rb @]. Much emphasis has been 
put on the attainment of the transition to a superfluid Cooper- 
paired state HGH 0. 

In a single-component Fermi gas at very low energies s- 
wave scattering is excluded by the statistics. Therefore, such 
a system behaves as a nearly ideal one. However, in this situ- 
ation other types of interactions, usually negligible, come into 
play. A good candidate is a dipole-dipole force. Its effects 
might be seen in the behavior of atoms with large permanent 
magnetic moments (such as chromium Jl2l[l3ll ) and should be 
evident in the case of heteronuclear molecules possessing per- 
manent electric dipoles (an example is ND3, which has been 
cooled in a Stark decelerator down to the millikelvin range 
11411 1. Another option is to induce electric dipoles with the 
help of strong DC fields 1 15] or by optical admixing of large 
dipole moments of Rydberg states 1 16]. 

The ground state of a harmonically trapped polarized 
single-component dipolar Fermi gas in the normal phase has 
been investigated in |17]. An interesting possibility of me- 
chanical instability of the system has been investigated in de- 
tail. Using the Thomas-Fermi theory (together with the Dirac 
correction) the authors have shown that when the critical val- 
ues of the total particle number N or of the dipole moment d 
are reached, the system collapses (for an analogous discussion 
of a dipolar Bose condensate see 1 16]). It has been shown that 
this behavior is governed by a parameter N*e, where 

s = (com 3 /h 5 f d 2 . (1) 

In Eq.Q co is the transverse frequency of an axially symmetric 
harmonic trap (the trap axis being parallel to the dipole polar- 
ization) and m is the particle mass. The critical value of N*s 



depends on the trap aspect ratio (3 = to z /oj (to z is the axial fre- 
quency of the trap). A detailed stability diagram is presented 
in Fig.l of Ref.|17]. Here we summarize its main features 
pertinent to the following discussion. The most striking find- 
ing is the existence of a characteristic value of the aspect ra- 
tio p 5.2 above which (i.e., in sufficiently oblate traps) the 
system is always stable, for any number of particles and any 
dipole moment. Away from this region the critical values of 
N^e, above which the system is unstable, tend to lie between 
2 and 10. Generally, the critical N*s is smaller for traps elon- 
gated in the direction of polarization, i.e., for /? < 1 (with an 
exception region for extremely stiff traps with /3 <§c 1, where 
the critical value of N*e grows again). As regards the exper- 
imental prospects of observing the predicted instability, one 
can say that it will definitely not affect weak magnetic dipoles 
I12U13I1 . However, polar molecules and strong electrically in- 
duced atomic dipoles may very well be located in the vicinity 
of the unstable region. 

At this point, a discussion of p-wave induced collapse in a 
single-component Fermi gas should also be mentioned 
Another interesting aspect of the physics of trapped dipo- 
lar fermions, an issue of superfluid pairing induced by the 
dipole-dipole interactions, has also been addressed recently 
inRefs.Gllii. 

The issue of collective excitations of degenerate Fermi 
gases is important and its studies have the same motivations 
as in the case of Bose condensates - the spectrum of collec- 
tive excitations determines the system dynamics in the low- 
energy regime. This way, the analysis of the excitations brings 
yet another important piece of information about quantum de- 
generate fermionic gases. A few works have been concerned 
with the problem. Amoruso et al. 12111 as well as Bruun 
and Clark 1 22] have investigated hydrodynamic excitations in 
non-interacting Fermi gases, mainly in spherically symmetric 
traps. Their work has been extended to the case of anisotropic 
trapping potentials by Csordas and Graham I23I1 . Some exper- 
imental work should be mentioned in the context of collective 
excitations - DeMarco and Jin have investigated spin excita- 
tions in a fermionic gas of atoms I24I1 . 

In the case of degenerate bosons, the spectrum of the collec- 
tive excitations strongly depends on the particle interactions. 
In the mean-field picture, this amounts to the dependence on 
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both the number of particles and the scattering length. As 
single-component degenerate fermions at low energies are vir- 
tually non-interacting, the problem of interactions is less im- 
portant and appears in two contexts only: either one deals with 
a two-component gas (where inter-component scattering takes 
place) 1 25] or, in the single-component case, one takes into 
account dipole-dipole forces. In this paper we study the latter 
case. 

In the hydrodynamic regime, i.e., when collisions ensure 
local equilibrium, the dynamics of a Fermi gas is governed by 
the hydrodynamic equations for the fields of velocity v(r, f) 
and density n(r, f): 
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where Vdd(r) is the dipole-dipole potential 



Vdd(r) = -=■ - 3 
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In the second of Eqs.(|3 one can recognize various com- 
ponents of the Thomas-Fermi formulation for trapped dipo- 
lar fermions 11711 : the first term on the right-hand side is the 
quantum pressure (microscopic kinetic energy) resulting from 
Fermi statistics; the next term describes the macroscopic part 
of the kinetic energy; the first term in the last line is the har- 
monic trapping potential and the last term in this line is the 
dipole-dipole interaction energy. 

In order to study low-lying collective excitations, in the fol- 
lowing we propose to merge the variational approach (used to 
find the ground state of a dipolar Fermi gas in Ref.[ 17]) with 
the time-dependent variational method used to study bosonic 
oscillations (see, e.g., Refs. l26l WA Eilo . A related ap- 
proach has been developed by Zaremba and Tso ll29ll to model 
a parabolically confined electron gas. 

First, we introduce the Lagrangian density 



mt) = - M r,0^ - K»(r, f )f 3 



dt m 20tt 2 
-|^»(r, i)(Vx(r, t)) 2 - l -mu> 2 [x 2 + y 2 + (J3z) 2 ] »(r, f) 

~J dr'n(r,f)Vdd(r-r')n(r',f). 
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where v(r, f) = (ti/2m) lim^^Vi - ^^i,^, and n(r, f) 
is the spatial one-particle density, x is usually termed the po- 
tential of the velocity field. The Lagrangian density X(r, f) 
is chosen such that the Lagrange equations for the conju- 
gate variables n and x> following from the Lagrangian L = 



J dr -C(r, t), reproduce the hydrodynamic equations (|2ji for the 
Fermi gas. 

From this point we proceed in a close analogy to the proce- 
dure employed for bosons in Ref.(28|- We introduce a Gaus- 
sian variational ansatz for the density n(r, t) : 



n(x,y,z, t) 
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where w v (t) are time-dependent widths, and a variational form 
for the phase: 
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where a n (t) and j3 n {t) are variational parameters - see also 
lEil - The linear part of the phase ansatz describes the 
drift of the Fermi fluid (a constant part of the velocity field) 
whereas the quadratic term corresponds to oscillations (a lin- 
ear part of the velocity). This kind of ansatz has proved to 
work very well in reproducing the static properties of trapped 
dipolar fermions 1 17]. It was also used successfully to model 
the dynamical behavior of trapped dipolar Bose-Einstein con- 
densates J27I E3l . This way, the problem is reduced to the 
study of a few Lagrange equations for the variational parame- 
ters w n (t), 770(f), a n (f), and/3, ; (f): 
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where c = 6 5/3 V375V /3 /25 is a numerical factor and 

I = ~J drdr'n(r,t)V M (r-r')n(r',t). (11) 

Again, we can recognize various terms of the Thomas-Fermi 
theory in Eas. ( 171101 : the first term of the right-hand side of 
Eq.Q corresponds to the quantum pressure, the second term 
comes from the trapping potential, and the last term results 
from the dipole-dipole interaction energy. 

As we see from Eq.©, the center-of-mass motion (de- 
scribed by 770(f)) is decoupled and follows simple harmonic 
oscillations, just as in the case of bosons 1 26, 28]. The only 
independent variables are the widths w n (t). Their dynamics, 
governed by Eq.Q, can be viewed as a motion of a fictitious 
particle in the following potential: 
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aspect ratio p 


P > 0.42 


P < 0.42 


lowest mode 


quadrupole (unstable for p < 5.2) 


N s e below the second crossing 


quadrupole 


e above the second crossing 


breathing (unstable) 



TABLE I: Lowest and unstable modes for harmonically trapped dipolar fermions. 



By diagonalizing the matrix of second derivatives of this 
potential (evaluated at the stationary point, i.e., for the ground 
state values of the widths - as calculated in Ref.| 17]) one can 
obtain the frequencies and the shapes of three low-lying eigen- 
modes of the cloud. These modes have the same structure as 
in the bosonic case - see Fig^ 



mode 1 



mode 2 



mode 3 



FIG. 1 : Graphical representation of oscillation modes given by the 
variational analysis. We refer to mode 2 as the breathing one and 
to mode 3 as the quadrupole one. Mode 1 will be called the 2D 
quadrupole mode. 



An interesting question in the analysis of the instability in- 
duced by dipole-dipole interactions concerns the identification 
of the mode that becomes unstable (which is indicated by its 
frequency approaching zero) at the critical value of the dipolar 
parameter N^s. We have analyzed this problem for different 
cylindrically symmetric traps. 

Let us first concentrate on cigar-shaped traps (J3 < 1). For 
0.42 < /J < 1, the lowest mode is always the quadrupole mode 
(mode 3 in Fig[3 and its frequency goes to zero as the crit- 
ical value of N^s is approached. At the characteristic value 
of ft - 0.42 two avoided crossings are present - see Fig|2] 
Away from the instability (i.e., for N^s < 2.40) and below 
both crossings the breathing mode (mode 2) has the highest 
frequency. The quadrupole mode is the lowest one and mode 
1 (2D quadrupole) is the intermediate one. At the first cross- 
ing, the nature of two higher modes is exchanged - now mode 
1 (2D quadrupole) is the highest one. This first crossing is 
present in all traps with (3 < 1. At the second crossing the 
breathing mode 2 acquires the lowest frequency and becomes 
unstable. For N^e slightly below the second crossing, the os- 
cillation of the quadrupole (lowest) mode is almost purely ax- 
ial (the eigenmode is (-6, +6, ~ 1), where 6 indicates a very 
small value). Just above the crossing, the shape of the lowest 
(breathing) mode remains very axial with no phase difference 
between the oscillations in the X and Y directions (the eigen- 
mode is now (+6, +6, 1)). This is always the case for/3 < 0.42. 



For smaller values of the trap aspect ratio, the positions of the 
two crossings merge (e.g., for/3 = 0.031 the crossings are at 
N*£ = 2.99, very close to the instability). For pancake traps 
(J3 > 1), the mode that drives the instability (which takes place 
as long as/? < 5.2) is always the quadrupole mode. For clarity, 
the results are summarized in TableU 

The time-dependent variational method employed thus far 
provides approximate solutions of the hydrodynamic equa- 
tions of the system. To confirm their validity, we have 
performed exact numerical simulations of the hydrodynamic 
equations 0. First, we find the ground state of the system by 
the propagation in imaginary time and use it as an initial con- 
dition for the simulations. Then, we perturb the trapping po- 
tential for a fraction of a trap period and let the system evolve 
while monitoring the resultant oscillations. Finally, we extract 
the frequency of the oscillations and identify the correspond- 
ing modes. The results are presented in Fig 13 and agree very 
well with the variational solutions. 




FIG. 2: Excitation frequencies of low-lying modes for the dipo- 
lar Fermi gas in a cigar-shaped trap with p = 0.42. The (unstable) 
mode 2 (breathing) is depicted with a solid line, while mode 1 (2D 
quadrupole) with a dotted line. Mode 3 (quadrupole) is marked by 
a dashed line. Lines depict the results of the variational analysis, 
while squares come from exact numerical solutions of the hydrody- 
namic equations. Inset shows a magnification of the second avoided 
crossing (mode 3 denoted by triangles and mode 2 by circles). 



To summarize, we have investigated the signatures of the 
mechanical instability in a single-component Fermi gas of po- 
larized dipolar particles confined in a harmonic trap by study- 
ing the spectrum of its low-lying hydrodynamic collective ex- 
citations. We have employed a time-dependent variational 
method, not previously used in the context of excitations of 
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dilute atomic Fermi gases, as well as exact solutions of the 
hydrodynamic equations of the system. We have identified 
modes driving the instability for various geometries of the 
trapping potential. In a rough classification, for oblate and 
slightly prolate traps (/} > 0.42) it is the quadrupole mode 
that drives the instability, while for cigar-shaped traps with 
/3 < 0.42 the breathing mode is an unstable one. Finally, 
one should note that the presented analysis describes the zero- 
temperature situation. It would be very interesting, although 
very challenging, to include the effects of temperature in the 
treatment. 
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